Correlated Phenotypic Transitions to Competence in Bacterial Colonies 
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Genetic competence is a phenotypic state of a bacterial cell in which it is capable of import- 
ing DNA, presumably to hasten its exploration of alternate genes in its quest for survival under 
stress. Recently, it was proposed that this transition is uncorrelated among different cells in the 
colony. Motivated by several discovered signaling mechanisms which create colony-level responses, 
we present a model for the influence of quorum-sensing signals on a colony of B. Subtilis cells during 
the transition to genetic competence. Coupling to the external signal creates an effective inhibitory 
mechanism, which results in anti-correlation between the cycles of adjacent cells. We show that this 
scenario is consistent with the specific experimental measurement, which fails to detect some un- 
derlying collective signaling mechanisms. Rather, we suggest other parameters that should be used 
to verify the role of a quorum-sensing signal. We also study the conditions under which phenotypic 
spatial patterns may emerge. 



Genetic competence is a cellular differentiation pro- 
cess, in which bacterial cells (such as those of Bacillus 
Subtilis) synthesize a specific system of proteins for the 
binding and uptake of DNA \l\ . Internalizing exogenous 
DNA leads to genetic transformations which presumably 
increases survival probability [2] . Under certain environ- 
mental conditions, a population of genetically identical 
cells (grown in homogeneous and identical environments) 
will at a given time be composed of two distinct subpop- 
ulations, a competent minority versus a non-competent 
majority, each with distinguished features. This has of- 
ten been attributed to bistability |31 IH [S] in the ge- 
netic regulatory network responsible for the competent 
state and indeed the necessary positive feedback has 
been found in the dynamics of one of the key regulatory 
molecules (comK) [BJ|7]- 

More recently, Elowitz and co-workers [7] have pro- 
posed that the competence transition is reversible and 
that a competent cell will return to its vegetative state 
after a characteristic time. This means that one should 
think of the network as being excitable (9j, with stochas- 
ticity (perhaps due to small molecular numbers) occa- 
sionally driving quiescent cells into the excited compe- 
tent state. This result is very interesting, as it represents 
one of the first examples where the inherent stochastic- 
ity of genetic networks [TUl ITTl IT^ is being put to 
functional use. One consequence of this new picture is 
that the competent subpopulation is not fixed but rather 
is continually fluctuating. Surprisingly, they have also 
claimed that this transition is made on the single-cell 
level, and is uncorrelated among different cells. This is 
supported by the measurement of the competence excur- 
sion times of neighboring cells. But in this work we show 
that this may not be the right parameter to measure, if 
quorum sensing is to be studied. The purpose of this 
paper is to present a putative quorum-sensing mecha- 
nism that induces short-ranged spatial anti-correlations 
between cells, and to suggest the appropriate way to ex- 
perimentally detect such a mechanism. This effect can 
be mediated by extracellular signaling molecules that are 



needed for the transition. In this scenario, a cell becom- 
ing competent will bias its neighbors against making the 
same choice. Yet, excursion times remain uncorrelated, 
as obtained experimentally. Our work is motivated by 
the many recently discovered instances in which signal- 
ing between bacterial cells fT4] are used to create colony- 
level collective responses to external environmental chal- 
lenges [inillS]. For instance, it is clear that the sporula- 
tion response in bacillus is spatially coordinated [T7]. At 
the end, we will discuss experimental tests of our predic- 
tions. 

We start by considering the molecular interactions un- 
derlying the competence circuit. The onset of compe- 
tence is stimulated by cell crowding and nutrient limita- 
tion. Detection of cell crowding is achieved by cell-cell 
signaling; various molecules are secreted and sensed by 
the cells, to monitor population density (quorum sens- 
ing) as well as other environmental conditions. In B. 
Subtilis, CSF (competence and sporulation factor) is a 
diffusible peptide, derived from the PhrC pheromone, 
which is secreted by the cells during growth. As its ex- 
tracellular concentration increases due to cell crowding, 
CSF is transported back into the cell TP, 120] where it 
influences a cascade of several genes, leading eventually 
to either sporulation or competence, depending on its 
concentration. For low concentrations (but above the 
transport threshold) it (indirectly) enhances the phos- 
phorylation of ComA, and consequently the expression of 
ComS. ComS competes with ComK for degradation by 
the MecA-ClpC-ClpP complex, thereby causing an indi- 
rect activation of ComK by ComS. ComK activates the 
transcription of several genes that encode DNA trans- 
port proteins and is considered as the indicator for the 
competent state. ComK activates its own expression and 
indirectly inhibits the expression of ComS. This double 
positive feedback loop can, as already mentioned, lead to 
a bistable response, i.e. two possible states with different 
gene expression. The negative feedback loop governs the 
escape from the competent state, as ComK degradation 
increases when ComS production is inhibited. 
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FIG. 1: Schematic representation of the competence circuit 
in B. subtihs. The autoregulatory positive feedback loop of 
ComK and the indirect repression of ComS are presented in 
the right-hand side. ComK and ComS compete for degra- 
dation by the complex MecA-ClpC-ClpP. The quorum sens- 
ing mechanism is presented in the left-hand side. Phrc is 
secreted by the cells and diffuses through the cell membrane. 
As its intracellular concentration increases, it diffuses back 
into the cell, initiating the competence pathway by regulat- 
ing the ComS production. ComK inhibits the production of 
SpoOA, which results in the inhibition of the cell crowding 
signal Phrc. Solid lines represent direct interaction, dashed 
lines represent indirect interaction. 



In ref. [71 [8] a model was proposed for the part of the 
circuit downstream from the external signal transducer 
ComA. Specifically, they proposed 
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where K and S represent the concentrations of the ComK 
and ComS proteins, (+ bj^) represent the minimal 
(fully activated) rates of ComK production, and Kq is 
the concentration of ComK required for 50% activation. 
The cooperativities of ComK auto-activation and ComS 
repression are parameterized by the Hill coefficients n=2 
and p=5. Similarly, the expression of ComS has maxi- 
mum rate 5s and is half-maximal when K = Ki . Degra- 
dation by the MecA complex affects both ComK and 
ComS; the nonlinear degradation term corresponds to 
such a competitive mechanism. Random fiuctuations in 
ComS expression are represented by a noise term rj(t). 
Degradation rates of K and S are normalized to unity 
and are assumed to be identical. The parameters can be 
chosen to give rise to an excitable system where there 
exists a stable fixed point plus two other intersections of 
the K-S nuUclines, a saddle point and an unstable fixed 
point. This case is presented in Fig. 2(a), where we also 
show a typical trajectory for which the noise has caused 
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FIG. 2: Phase planes diagram for the ComK-ComS system 
(Eqs. (l)-(2)). Null clines are presented for ComK and 
ComS, and the system fixed points, (a) The excitable system, 
obtained from Eqs. (l)-(2) with the parameters: a=0.004, 
6a'=0.08, &s=0.8, fco=0.2, fci=0.222. The noise RMS ampli- 
tude is 3.5 and At =0.01. The fixed points are (from left to 
right) a stable point, a saddle point and an unstable point. A 
characteristic trajectory is presented for a single cell, staring 
from the vegetative state (low ComK). (b) The monostable 
system, for Eq. (3) and g=0.9, 6^=0.6, Cs=0.2, Sq=1, Dq=1 
(other parameters as in (a)). The only fixed point refers to 
the vegetative state. 



a large excursion, i.e. has excited the competent state. 
Eventually the cell returns to the vegetative stable state, 
due to the repression of ComS by ComK; when the value 
of S decreases below a threshold level, K is more rapidly 
degraded and can no longer auto-activate its production. 
One can obtain similar results with more biologically re- 
alistic noise forms [H [8l [21] . 

In this model, each cell was taken to be independent 
of the other cells. This assumption cannot be completely 
correct, however, since the external enabling signal act- 
ing through ComA comes from the rest of the colony. 
Obviously, the functional utility of having a competent 
subpopulation would clearly be enhanced if this subpop- 
ulation were spatially dispersed. In order to consider the 
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possible spatial effect of quorum sensing, we introduce an 
extension of the model, and take into account a putative 
biological mechanism (see Fig. 1), in which the basal 
production level of ComS depends on a quorum sensing 
pheromone Q. The motivation for this particular model is 
a microarray analysis |22j which has shown that the gene 
SpoOA, which is upstream of the PhrC secretion mecha- 
nism, is down regulated when ComK is high. Thus, the 
secretion of CSF, which is a regulator for competence 
onset, is effectively inhibited by ComK. Specifically, we 
assume that 



Now K and S are functions of space as well as time, and 
Q depends on the cells' state in the following manner: 

^ = i tanh {A{K^ax ~K) + \) + DV^Q - SqQ (3) 

The production of Q is maximal as K —> , namely nor- 
mal growth, and it decreases as the cell becomes compe- 
tent. Kraax IS choscn to bc the threshold value of K, used 
to define the competent state. Q diffuses with diffusion 
coefficient D , and is degraded with a constant rate 5q. 
The parameters are chosen to describe a bistable system 
for Q = 1, an excitable system for Q close to but less 
than one and a monostable one if Q << 1. In fact, for 
the system described in Fig. 2(a), Q level as low as 0.9 
is sufficient to drive the system to the vegetative monos- 
table regime, as can be seen in Fig. 2(b). This is a key 
point, as very small changes in the external signal (e.g. 
the level of CSF) may drive the cell into a completely 
different regime. Therefore, even a 2-fold change in the 
production of SpoOA, and as a result in the secretion of 
CSF, may have a substantial effect on a possible compe- 
tence transition. 

Our goal in this model is a semi-quantitative evalua- 
tion of correlations in the competence status of nearby 
cells. We will study this question in one spatial dimen- 
sion (with periodic boundary conditions), with a cell at 
every one of the 5000 grid points. In each explicit time 
step, the values of K and S are updated for each cell 
and then Q is updated in the entire space. A typical 
simulation run of our model starting from all vegetative 
cells is shown in Fig. 3. There are clear anti-corelations 
between cells, due to the effective inhibition created by 
a reduction of Q emission by the competent cells. To get 
a more quantitative measure, we measured the spatial 
correlation function 

Cd{i) = (4) 

where we average over time, a is the standard deviation, 
and the subscript refers to a specific cell i. We then aver- 
age over all cells i to obtain an estimate Cd- In Fig. 4a we 
present the function Ci , i.e. the correlation between ad- 
jacent cells, versus the diffusion coefficient of the quorum 
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FIG. 3: Time development of a colony of cells. Color code 
represents the ComK level, from low ComK level (vegetative 
state, black) to high ComK level (competent state, white). 
The system parameters are as in Fig. 1(b). All cells are 
initially vegetative, and competence events occur with an av- 
erage waiting time of 45 hours. The average excursion time 
in the competent state is about 5 hours. Adjacent cells rarely 
become competent simultaneously, due to effective mutual in- 
hibition. 



sensing pheromone Q. For a small diffusivity, the inhibit- 
ing effect is negligible as the diffusion length is small com- 
pared to the inter-cell distance. As D is increased, the 
mutual inhibition is increased and the correlation is neg- 
ative, as expected. The maximal effect on nearest neigh- 
bors is when the characteristic length A = ^/ Dt, given by 
the diffusivity and the degradation time t = is com- 
parable to the inter-cellular distance. For the given pa- 
rameters, the maximal effect on adjacent cells, i.e. A=l, 
is obtained for D = \. For an inter-cellular distance of 
1/im, the highest mutual inhibition on nearest neighbors 
will be obtained for D — 1/t {^.m? /sec). 

In the experimental data of [7 , sister cells showed no 
correlation between their excursion times, once they be- 
come competent. This led to the assumption that there 
was no important spatial coupling. However, this conclu- 
sion might be inaccurate, as there might be correlations 
between the cycles of neighboring cells but not in their 
excursion times. In our simulations, excursion times only 
slightly changed with the diffusion coefficient (about 10% 
- data not shown). And, no correlation was observed be- 
tween the excursion times of neighboring cells, as can 
be seen in Fig. 4(b), which is very similar to the re- 
sults of [7 . However, there are notable correlations in 
the times of competence onset of neighboring cells, as 
demonstrated in Fig. 4(a). This occurs because the Q 
field has little effect on the competent state itself, but 
is required to be high for the noise-induced excitation. 
Therefore, the excursion times measurement is mislead- 
ing and does not definitively indicate whether there is, 
or there is not, a quorum sensing mechanism which influ- 
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FIG. 5: A striped colony of cells simulation results. The cell 
state is represented by color code (see Fig. 3). (a) For D=l, 
the stripes pattern is unstable and the cells return to their 
vegetative state, (b) For D=2, the stripe pattern is stable. 



FIG. 4: (a) Correlation between ComK values of adjacent 
cells as a function of diffusivity. The mutual inhibition, as 
identified by negative correlation, vanishes for small diffusiv- 
ities as well as for large diffusivities. (b) Cumulative his- 
tograms of differences between excursion times, of the entire 
colony (solid line) and adjacent cells (markers). 



ences the competence transition. Clearly, our predictions 
could be tested within the same experimental framework. 

The existence of diffusing inhibition in this system sug- 
gests that under certain circumstances the colony should 
be able to support a patterned Turing state consisting 
of periodically repeating regions of competent and non- 
competent cells. There are two distinct possibilities; the 
system might have a Turing instability which eliminates 
the vegetative state completely or it might be multi- 
stable, meaning that patterned states co-exist with the 
uniform one. We have found no evidence of instabil- 
ity, but multistability can exist for large enough diffusion 
constant. In Fig. 5 we consider the temporal evolution 
of a colony with an initial stripe pattern, in which ev- 
ery tenth cell is competent. The values of the relevant 
parameters are chosen to be: Ofe — 0.004, 6^ = 0.14, 
/co = 0.2, hs = 0.58, fci = 0.222 and = 0.1. For 
these values the K-S system of equations is bistable for 
Q = 1 and vegetative monostable for Q < 0.6. When the 



diffusion of Q is slow (Fig. 5(a)), the local decrease in 
Q shifts the local K-S system of each competent cell to 
the monostable, vegetative state, the cells revert to veg- 
etative and the pattern vanishes. But when diffusion is 
faster (Fig. 5(b)), the averaged value of Q can lie above 
the monostability threshold, and the pattern is stable. 

We can get a better understanding of the stripe pat- 
tern and of the lack of a Turing instability by explicitly 
evaluating the Q field. If we approximate the tanh in cq. 
3 as a step function, we find that in steady-state. 



Qix) 




(5) 



where the sum ranges over competent cells. This geomet- 
ric series can easily be calculated for any specific x; the 
result is shown in Fig. 6 as compared to the simulation 
data. We see that the pattern is possible only because the 
self-inhibition is weakened sufficiently by the diffusion so 
as to allow the competent sites to remain in the bistable 
range. Paradoxically, the competent cells are the most in- 
hibited, which is presumably why this pattern would not 
spontaneously appear from the uniform vegetative state. 
Studying this experimentally would necessitate prepar- 
ing the system with an externally imposed pattern of 
CSF and then determining whether the pattern remains 
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behavior of a colony of B. Subtilis cells. The condition of 
each cell is described by the concentrations of two pro- 
teins, ComK and ComS. The dynamics is described by a 
set of two rate equations, with a double-feedback loop. 
The production rates effectively depend on the cell den- 
sity in the entire colony, sensed by a diffusible pheromone. 
This pheromone is likely to depend on the protein ComK. 
This additional mechanism creates an effective inhibition 
between neighboring cells, leading to anti-correlation in 
the cells competent/ vegetative cycles. Excursion times of 
neighboring cells arc found to be uncorrelated, due to the 
self-excitatory mechanism of ComK, also in agreement 
with experimental results. These findings can be exper- 
imentally verified by the direct measurement of ComK 
concentration in the cells, as was done in ref. [7] via flu- 
orescence techniques. 



FIG. 6; The steady state profile of Q(x), analytic solu- 
tion(solid line) and simulation results (asterisks) 



in place even after the external pattern is removed. 
To summarize, we present a new model for collective 
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